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ABSTRACT 

If the dark matter particle is a neutralino then the first structures to form are cuspy 
cold dark matter (CDM) haloes collapsing after redshifts z « 100 in the mass range 
10~ 6 — 10~ 3 M Q . We carry out a detailed study of the survival of these micro-haloes 
in the Galaxy as they experience tidal encounters with stars, molecular clouds, and 
other dark matter substructures. We test the validity of analytic impulsive heating 
calculations using high resolution A-body simulations. A major limitation of analytic 
estimates is that mean energy inputs are compared to mean binding energies, instead 
of the actual mass lost from the system. This energy criterion leads to an overestimate 
of the stripped mass and underestimate of the disruption timescale since CDM haloes 
are strongly bound in their inner parts. We show that a significant fraction of material 
from CDM micro-haloes can be unbound by encounters with Galactic substructure and 
stars, however the cuspy central regions remain relatively intact. Furthermore, the 
micro-haloes near the solar radius are those which collapse significantly earlier than 
average and will suffer very little mass loss. Thus we expect a fraction of surviving 
bound micro- haloes, a smooth component with narrow features in phase space, which 
may be uncovered by direct detection experiments, as well as numerous surviving cuspy 
cores with proper motions of arc-minutes per year, which can be detected indirectly 
via their annihilation into gamma-rays. 

Key words: cosmology: theory - dark matter - galaxies: formation - gamma-rays: 
theory - methods: numerical 



1 INTRODUCTION 

If dark matter is composed mainly of the lightest super- 
symmetric partner particle, the neutralino, the first self- 
gravitating structures in the Universe are Earth-mass haloes 
forming at high redshifts (|H qfmann . Schwarz fc Stocked 
120011 iDiemand. Moore fc Stadell l2005lh As many as 10 15 
could be within our Galactic halo today. These abundant 
cold dark matter (CDM) micro-haloes have cuspy density 
profiles that can withstand the Galactic tidal field at the 
solar radius. The numbers of such haloes that lie within the 
vicinity of the solar system depends on how many survive the 
complex merging history of early hierarchical structure for- 
mation. iV-body simulations of CDM satellites indicate that 
tightly bound cusps are very stable against tidal stripping 
(Kazantzidis et al. 2004), and therefore dense micro-haloes 
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accreting late onto more massive structures may survive rel- 
atively intact. The exact distribution of dark matter in the 
solar vicinity is important for direct and indirect dark mat- 
ter detection experiments. 

Substructures that survive the merging process will 
experience continuous perturbative encounters with stars, 
molecular cl ouds, and other dark ma t ter su bhaloes. As 
discussed in iDiemand. Moore fe Stadel (I2005D . we expect 
that these encounters lead to some mass loss but that 
the cu sps of most micro-h al oes remain intact. Recent stud- 
ies bv | Z hao et alJ i2005 allbl). iGre en fc G oodwin! (2006) and 

have raised the 



iBerezinskv. Dokuchaev fc Eroshenkd i200 
question whether these first haloes would be completely dis- 
rupted by close encounters with stars. Crossing the Galactic 
disc w ould also cause additional tidal heating. iMoore et alJ 
(2005) argued that the analytical impulse approximation 
and the semi-analytic models used in these studies may not 
fully describe the disruption of the high-density inner cores. 
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Particle orbits deep in the cusp may remain adiabatically 
invariant to the perturbations and preserve the structure 
of the cusp. Only direct numerical simulations can describe 
these complex dynamical processes. In this paper, we use 
several sets of high resolution iV-body simulations to test 
the validity of analytical heating models. 

An important factor in the survival statistics of micro- 
haloes is how many survive sim ilar-mass mergers during the 
build -up of the Galactic halo jPiemand. Kuhlen fc Madaul 
2006). Even if only a few percent survive the hierarchical 
growth, many micro-haloes would still lie within one parsec 
from the Sun. Their dense cuspy cores would be sources of 
gamma-ray emission due to self-annihilation, which could be 
uniquely distinguished by their high proper motions on the 
sky of the order arc-minutes per year (Moore et al. 2005, 
Koushiappas 2006). 



2 HEATING BY STARS IN THE SOLAR 
NEIGHBOURHOOD 

There are various ways to define a virialized halo. The ap- 
proach often used in cosmological simulations, which we 
adopt here, is that dark haloes virialize when their aver- 
age density equals A = 200 times the mean density of the 
Universe, p(z) = 3f2 #o (l+z) 3 /8ttG. Here H is the Hubble 
constant, fio is the matter density parameter, and z is the 
redshift of virialization. Then the virial radius of the halo, 
defined by the relation M V1I = ^-R 3 {l Ap(z), is 



i?vir = 0.31 (l + z)" 



M vir 



10" 6 M 



1/3 



pc. 



(1) 



for Q = 0.3 and H = 70 km s" 1 Mpc" 1 . The virial velocity 
is defined by the relation V v 2 lr = GM v i r /7? v i r : 



Kir = 12 (1 + z 



,1/2 



Mvir 



1O" 6 M 



1/3 



(2) 



These parameters determine the binding energy of the 
haloes, which can be expressed usin g the half-ma ss radius 
of the system: E b w 0.2GM vir /R 1/2 iSpitzerlll987b . Density 
profiles of dark matter haloes in cosmological simulations 
are often described by the NFW model with a concentra- 
tion parameter, c. For c < 10, the radius containing half of 
the virial mass is approximately R1/2 ~ (5c) _1//4 i? v ir- High- 
redshift haloes have typically low concentrations, such that 
R1/2 ~ 0.5-Rvir- Therefore, the binding energy of first haloes 
is E h w QAV^- 

As these small haloes merge into larger systems, two 
effects may modify their structure: tidal truncation by the 
host galaxy and tidal heating by massive, fast-moving per- 
turbers (stars, molecular clouds, other dark matter substruc- 
tures). 

In the vicinity of the Sun, the matter density is 
dominated by stars, which we assume to have the same 
mass m* = 0.7 Mp>. The stellar mass d ensity is m*n» w 
O.IMqPc" 3 l|Binnev fc Merrifieldl Il99sl. . which is a half 
of th e total density of the disc calculated from the Oort 
limit flBahcalll llbs^ . In order to remain gravitationally self- 
bound, the micro-haloes must have an average density above 
roughly 2m,n 4 . 



Fast encounters with massive perturbers increase the 
velocity dispersion of dark matter particles and reduce a 
halo's binding energy. A distant encounter at an impact pa- 
rameter 6 with a relative velocity V vc \ increases the energy 
per unit mass on the average by 

1 / 2Gm, * 2 



ABi(6) 



2 U 2 Kci 



(3) 



where (V 2 ^ ~ Rf/2 is the ensemble average of the particle 
distance squared from the centre of the micro-halo. 

At very small impact parameters, 6 < bi, a single en- 
counter would be sufficiently strong to unbind the whole 
halo: A_Ei(6i) = Eb- As we show later in section 2] a small 
central part always survives even such a strong perturba- 
tion, apart from direct collisions with 6 = 0. Nevertheless, it 
is instructive to define the disruptive encounter threshold, 
which is given by 

Gm,R vi ~\^ 2 



61 



/ Gm, R vir V 

3C V KclKir ) 



0.2 (l + z) 



-3/4 



pc, 



(1) 



KelKir 

where a c w 0.96 (c/3) _1/ ^ 8 . Equation 10 is strictly valid only 
in the tidal approximation, b 3> -Rvir- An encounter at 61 
falls in that regime for redshifts z < 50, which is appropriate 
for our consideration of the micro-haloes. 

The number of encounters over time t as a function of 
impact parameter is dN cnc (b) — n t V rc it2nbdb, where n* is 
the number density of stars. We can obtain the cumulative 
effect of multiple non-disruptive encounters by integrating 
over the impact parameter: 

AStid = f 

Jb 1 



AE^^^db 
db 



4 G 2 m 2 ,Rl ir n*t ( 1 



V- 



ti 2 b 2 



(5) 



The upper limit of integration is set by the condition that 
the encounter is impulsive, i.e. the duration of the encounter 
6/Vrei is shorter than the orbital time of particles in the halo, 
Rvii/V v ir. The maximum impact parameter is given by 



/ bmax \ 



al ffi ^T ir » 1. 



(6) 



bi J c Gm»V vi 
The ratio of the tidal heating energy in non-disruptive en- 
counters to the binding energy is 

AE ti d 2 Gm t n t tR viT 

— — =a c 7r . (7) 

Hb Vvir 

We can also calculate the effect of disruptive encounters, 
with b < 61. The number of such encounters is simply 

, , > 1 2 t r 1 2 Gm*n*tR v ir , . 

N cnc {< bi) = Tvb^tVrcit = a c -K — . (8) 

V vir 

This number is the same as equation Q meaning that the 
cumulative effect of disruptive encounters is the same as that 
of non-disruptive encounters. The total disruption probabil- 
ity, iVtot, is then twice that given by equation J7[. 

To calculate this disruption probability, we note that 
while stars in the solar neighbourhood move on approxi- 
mately circular orbits around the Galactic centre, small dark 
matter haloes would be moving on isotropic orbits inclined 
with respect to the Galactic disc. Their expected vertical 
velocity is V z ~ 200 km s _1 . The crossing time of the disc 
with a scale height of H = 0.2 kpc is 2H/V Z = 2 x 10 6 yr. In 
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the solar neighbourhood, haloes would cross the disc every 
10 8 yr and have about 100 crossings in the Hubble time. The 
total amount of time the haloes would spend in the region 
of high stellar density m,n, is then td ~ 2 x 10 8 yr. The 
total disruption probability is 

jV tot = 2iV cnc (<6 1 ) 

= (l + zy 3 '*( m*n, W t d \ 

V 131 >/ \0.1M e pc-3 J \2 x 10 8 yr/ ' V ; 

Therefore, the haloes virialized after redshift z — 130 should 
suffer significant mass loss by passing stars in the solar 
neighbourhood. Due to biased halo formation typical sub- 
haloes in the solar neighbourhood come form 2a fluctuations 
(Diemand, Madau & Moore 2005), i.e. they virialize at half 
the expansion factor (or twice the (z + 1) value) than typical 
haloes of the same mass in the field (i.e. la peaks). A forma- 
tion time of z = 130 corresponds to about a 3a peak. Such 
early structure formation is not uncommon in dense environ- 
ments, for example the small, over dense region simulated 
in Diemand, Kuhlen & Madau (2006) already contains 845 
micro-haloes at z = 130. A fraction of about 20% of the 
local mass comes from peaks above 3a (Diemand, Madau 
& Moore 2005), implying that approximately 20% of the lo- 
cal subhalo population should therefore not suffer significant 
mass loss. 



where "j^ b is the number density of subhaloes /. Taking an 
NFW model for the smooth component of the Galactic halo 
and restricting our analysis to the inner part of the halo, 
r < r s w 20 kpc, we find the subhalo's truncation radius 
fi sll b ~ Vsifr/rs) 1 ' 3 . The density of subhaloes, assuming 
they have not been completely disrupted, is 

rin sub ^ IP- 2 /" 2 
df ~ 47n ? (c)r 3 V { ' 

where g(c) = In (1 + c) — c/(l + c) ~ 1.6 for a concentr ation 
parameter c = 12 jKlvpin. Zhao fc Somervillej 12002). For 
the Galaxy, a Hubble time corresponds to t ~ 5-Rhost /Vhost ■ 
Thus for the inner halo, 

Ag tid (/) _ 5 x 10- 2 c 5 / tV host \ i^Xost 
E b 4g(c) \5R hoB J Rl oat V v l 

x r 2/3 (^" 5/3 - (13) 

Integrating over all subhaloes, / < 0.01, we find 

*f-« 0.063 (LY 5/3 . (14) 
&b V r s J 

Thus, mini-haloes may be disrupted by repeated encounters 
with more massive haloes within r < 0.2r s m 4 kpc from the 
centre of the Galaxy. 



3 HEATING BY DARK MATTER 
SUBSTRUCTURE 

Virialized, self-gravitating subhaloes within larger haloes 
(the substructure) will also kinematically heat and disrupt 
their sm all cousins |c.f. Boily et al. (20 04)1. iV-body simula- 
tions dDiemand. Moore fc Sta dcl 2004) show that the num- 
ber of subhaloes with masses above M within a host of mass 
Afhost scales as 

iV 8ub (>M/M w) «( IF 4-)". (10) 

Since stars occupy only a small fraction of the volume of 
their host haloes, it is important to consider if the tidal 
heating by dark matter subhaloes can disrupt a significant 
fraction of micro-haloes. 

The analysis of section [5] can be generalised for per- 
turbers with a range of masses, M V1I < M < 10~ 2 Mhost- 
Let / = M/Mhost be the dimensionless subhalo mass. The 
threshold impact parameter, at which a single encounter 
with subhalo / is disruptive, is b\(f) » /VhostRjiostR^vir/Vvh-, 
where we assumed the relative velocity to be the virial ve- 
locity of the host halo, I^-el ~ Vhost- However, for most 
subhaloes this impact parameter is smaller than their size, 
Rsuh ~ '"[M su b/3A/host('")] 1 ' /3 , which is determined by tidal 
truncation at distance r from the centre of the host halo. 
Tidal approximation applies only at b > 6 m i n = i? S ub ■ There- 
fore, most encounters will be non-disruptive. 

The cumulative heating by multiple non-disruptive en- 
counters with subhaloes of mass M S ub = fMhost is [see equa- 
tion JjjJ]: 

AgtidCf) _ 7vG 2 Ml b R 2 vir t dn suh 

E b V hoBt V v lbl in df ' ( > 



4 NUMERICAL TESTS OF THE IMPULSIVE 
APPROXIMATION 

In this section we test the response of a CDM micro-halo to 
repeated impulsive encounters using N-body calculations in 
order to test the validity of the impulse approximation, and 
to study in detail how the internal structure of the micro- 
haloes evolve with time. 

The initial state for the micro-halo is an equilibrium 
profile with the same structural parameters as found by 
iDiemand. Moore fc Stadell (120051) at z = 26, the epoch at 
which such structures are typically accreted into larger mass 
systems. T his halo obeys a c uspy density profile, the general 
a, (3, 7 law teernauisd(l99Ch : 

P(r) = P - — (r ^ i? mi cro), (15) 

(^) 7 [l + (7r)1^ 

with a = 1.0, P = 3.0 and 7 = 1.2. The mass of the 
halo is M m icro = 10~ 6 Mq within the z = 26 virial ra- 
dius -Rmicro = 0.01 pc. The concentration parameter is low, 
Rmicro/r s = 1.6, typical of micro-haloes in the field at z = 
26. Some experiments we repeated with micro-haloes hav- 
ing concentrations of 3.2. The typical local subhalo forms 
earlier (by about a factor of two in redshift, see Diemand, 
Madau & Moore 2005) than the average micro-halo in the 
field. Therefore the typical local subhalo might be twice 
as concentrated and more robust against mass loss. How- 
ever, to be conservative we use the low concentration of 1.6 
throughout this paper unless stated otherwise. For numeri- 
cal stability of the profile, we make a realization of this halo 
ex tending to approximately 4 fl m i C r o usin g the techniques 
of iKazantzidis. Maeorrian fc Moord i2004f) . At r > -Rmicro, 
the density profile falls off exponentially as exp(— r/r^ ccajy ), 
with rdecay = 0.3 -Rmicro- The total mass of the system is 
therefore 1.81 M m i cro . We use 10 6 particles of equal mass, 
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m p = 1.81 x 10 _1 Mq. The force calculations have a soft- 
ening length of 0.005 -Rmicro- 

We then subject the equilibrium micro-halo to a se- 
ries of impulsive encounters with a star of mass m* = 
0.7 Mq, the mean mass per star in the disc of the Galaxy. 
First we run six simulations, which differ in the mini- 
mal distance from the star to the centre of the micro- 
halo. For these six simulations the impact parameters are 
b = 0.005, 0.01, 0.02, 0.05, 0.1, 1 pc. In all runs the star moves 
with the relative velocity V ro i = 300 km s^ 1 . The initial sep- 
aration of the star and the halo along the direction of mo- 
tion is three times the impact parameter or three times the 
micro-halo radius of the halo, whichever is the greater. After 
the star reaches the point of closest approach, we let it move 
away the same distance from the halo. Then we remove the 
star and let the system evolve in isolation for 3 x 10 s yr, 
which corresponds to 20 crossi ng times of the halo. Si milar 
experiments date back to e. g. lAeuilar fc White! (Il985h . 

Each encounter increases the internal energy of the 
micro-halo. Following the perturbation, the system under- 
goes a series of virial oscillations (contraction and expan- 
sion) until th e potential relaxes into a new equilibrium 
configuration jGnedin fc Ostrikedll99Eh . Depending on the 
strength of the perturbation, this potential relaxation takes 
between 10 and 20 crossing times of the halo. The final con- 
figuration has a lower binding energy and some particles 
may escape the system entirely. 

Fig- shows the energy change per particle for a flyby 
at b = 0.02 pc at two different times: directly following the 
encounter and after the potential relaxation. The duration 
of the encounter, r = 2b/V ve \ ~ 130 yr, is much shorter than 
the dynamical (crossing) time of the particles in the micro- 
halo, tdyn ~ 1.5 x 10 7 yr. Therefore, we expect that the tidal 
heating can be calculated in the impulsive approximation 
[c.f. equation JJ)]. For ensemble-average of stars with initial 
energy E, the energy per unit mass increases by the amount 

This prediction is plotted next to the numerical result 
(squares) in Fig. and agrees with it very well. 

Subsequent potential relaxation reduces the depth of 
the potential w ell of the system, lea d ing to another effective 
energy change. iGnedin fc Ostrikerl (I1999T) found that this 
additional energy change can be approximated as a constant 
fraction of the initial potential, $i: 

A£ pot (r) = c (-$;(*•)), (17) 

where the constant c is such that the sum of AE pot (r) over 
all particles is twice the initial energy change of the system, 
AEi(b), as required by the virial theorem. The final energy 
difference is {AE}+AE pot . This prediction is plotted next to 
the numerical result (triangles) in Fig.Qand again provides 
a very good fit. 

Fig. [5] shows the energy changes immediately follow- 
ing the encounters at different impact parameters, b. The 
analytical formula lib! provides a good description of the 
numerical results, except in cases of extremely strong per- 
turbations when the energy changes by more than 100%, 
(AE) > \E\. 

Fig. [3] shows the final energy redistribution after po- 
tential relaxation, for encounters with different impact pa- 




Figure 1. Energy change per particle as a function of radius im- 
mediately following the perturbation (squares) and after potential 
relaxation (triangles), for the encounter with b = 0.02 pc. 



rameters, b. Equations l|16^ and l|17^ describe the effect very 
accurately. 

The change of the micro-halo potential leads to the 
change of the density profile. Particles that gain enough en- 
ergy escape the system form unbound tidal tails. The final 
density profiles for the encounters with different impact pa- 
rameters are shown in Fig. [I] 

The amount of mass stripped from the halo depends on 
the definition of the bound mass. The density of the outer 
halo profile extends as r -3 beyond the nominal micro-halo 
radius, and all of the particles are initially bound. We use 
two practical definitions, (i) We have defined an effective 
maximum radius, i? ma x = 4i? m i cro , beyond which all parti- 
cles are assumed to be lost of the micro-halo. In practice, 
this radius can be set by the external tidal field, (ii) We 
have also defined an effective tidal potential at that radius, 
$t = 3?(-Rmax)- We use this tidal potential to construct an- 
other definition of unbound particles, as those with E > $t- 
After the potential relaxation, the new potential <E> is used 
to define <E> t at the same fixed radius -R max - 

Fig-OOshows the change of the total energy of the system 
immediately following the encounter, for all particles within 
-Rmax- It is well described by equation JjjJ, which is plotted as 
a solid line. Since the density profile of the system continues 
beyond the nominal micro-halo radius, the average radius 
for all particles is (r 2X j = 1.63 -R^ icro . We used this value for 
the plotted prediction. 

Fig. HJ shows the mass loss as a function of impact pa- 
rameter, using the two definitions based on the position and 
energy criterion, respectively. For large impact parameters 
(weak perturbations), the position criterion indicates sys- 
tematically lower mass loss than the energy criterion. There- 
fore, some particles within 7? max may be unbound at the end 
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Figure 2. Energy change per particle as a function of radius 
immediately following the perturbation, for encounters with dif- 
ferent impact parameters b. 



Figure 3. Energy change per particle as a function of radius 
after potential relaxation, for encounters with different impact 
parameters b. 



of the simulation. In the strong perturbation regime, both 
criteria give similar results. 

While the total energy change of the system can be 
computed with sufficient accuracy using the impulsive ap- 
proximation, the amount of mass lost cannot. Using our nu- 
merical simulations, we seek to establish a practical relation 
between AM/M and AE/E h . We find that the following- 
equation provides a good fit to the numerical results: 

AM _ 
~M~ ~ 



2.1 



AE 



(18) 



For weak perturbations the mass loss scales as the energy 
change, AM/M ~ 0.5AE/ Eb- For very strong encounters, 
the mass loss asymptotically approaches unity. Note how- 
ever, that even for very small impact parameters, when 
AE/Eb 3> 1, a small fraction of the mass always remains 
bound, « 2(AE/E b y 1 . 

A note on notation. Strictly speaking, the tidal approxi- 
mation which was used to derive equation Q is only valid for 
b 3> -Rmicro • At smaller impact parameters the energy change 
does not scale as b~ 4 , bu t it can be calculated in the oppo- 
site asymptotic limit (e.g, ]Moore |l993l : ICarr fc Sakellariadoul 
Il999l : iGreen fc Goodwin! l2006fl . We take an alternative ap- 
proach by parametrising the mass loss in our numerical sim- 
ulations (equation 1181 1 based on the formal extrapolation 
of equation J3J to all values of b. 

As Fig. 2] shows, most of the mass remaining bound 
to the micro-halo after strong perturbations is concentrated 
near its centre. It is therefore interesting to calculate the 
fraction of lost mass that was initially contained within the 
power law density cusp, at r < r s . Fig. Q shows that this 
fraction scales as AM cusp /M cusp « (1 + 5(AE/E h )~ 1 )~ 1 . 

We also performed another set of simulations, by re- 
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Figure 4. Density profile of the micro-halo in a new equilibrium, 
after encounters with different impact parameters, b. The arrow 
indicates the micro-halo radius. 



peatedly perturbing the micro-halo with the same star, with 
the same relative velocity and at the same impact param- 
eter b = 0.02 pc. After the halo has relaxed following the 
first encounter, we move the centre of mass of the remaining 
halo to the origin of the coordinate system and let a star 
pass by in exactly the same way and again let it relax, and 
then repeat the encounter a total of ten times. 

The density profile of the micro-halo after n encounters 
is shown in Fig. |H| Every encounter heats the system and 
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Figure 5. Total energy change of all particles immediately fol- 
lowing the perturbation, as a function of impact parameter {filled 
squares) and the analytical prediction in the impulse approxima- 
tion {solid line). 




Figure 7. Mass lost inside the scale radius, r < r s , after the 
potential relaxation, as a function of impact parameter. 
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Figure 8. Density profile of the micro-halo after successive en- 
counters with a star with the same mass and orbital parameters 
{b = 0.02 pc), marked by the encounter number n. 



Figure 6. Mass loss of the halo as a function of the impact param- 
eter b, determined using either the energy criterion, E > $t {open 
squares: immediately following the perturbation, filled squares: af- 
ter the potential relaxation) or the position criterion, r > i? ma x, 
after the potential relaxation {triangles) . 



strips some mass, but each successive encounter is less and 
less effective. Fig. [5] shows the cumulative mass loss after 
each such encounter. It can be described by the following 
function: 

^(n) = l-exp(An fl ), (19) 

with A = -0.34 and B = 0.87 for c = 1.6 and A = -0.23 and 
B = 0.81 for c = 3.2, which are also shown in Fig. 

We fitted each of the eleven density profiles, which are 
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Figure 9. Cumulative mass loss of different micro-haloes after a 
number of identical encounters with b = 0.02 pc. 

shown in Fig. |H| as well as the corresponding eleven den- 
sity profiles for the halo with c = 3.2, which are not shown, 

with equation 1151 . To make the fits we used the Levenberg- 
j ■ i ° 

Marquardt method (Marquardt 1963) in log- log space keep- 
ing 7 fixed at 1.2. The variation of a,/3,r s and po from en- 
counter to encounter is as follows: a increases from 1.0 to 
1.5, /3 increases monotonically from 3.0 to 7.0, r s stays con- 
stant around 7.0 milli-pc for the c = 1.6 halo and around 4.0 
milli-pc for the c = 3.2 halo and finally po oscillates around 
20 GeV/cm 3 for c = 1.6 and decreases from 100 down to 
50 GeV/cm 3 in case of the other halo. The resulting profiles 
could now be used to calc ulate the net flux c oming from 
neutralino annihilation fe.g. iKoushi appas 200JJ) via: 

/oo 
Anr 2 p(rfdr (20) 

We have summed up the dependence of the flux on neu- 
tralino mass and interaction cross section in the constant 
k. The lower bound r m i n is defined as the central region of 
the micro halo, in which the neutralinos already annihilated 
each other. The required number density for this to happen 
can be estimated with the help of 

th = — (21) 
nav 

where th ~ 13 Gyrs is the Hubble time, av ~ 10 -30 cm 3 s _1 
is a typical cross section and n is the number density of 
neutra linos. For more details see lCalcaneo-Roldan fc Moord 
(1999). The minimum radius can now be computed from 
comparing this minimum number density with the density 
profile in Fig.|H| Assuming a neutralino mass of 100 GeV and 
deploying the above mentioned density profile, r m i n comes 
out to be 1.6 x 10 -14 pc. Fig. 1101 shows then the resulting 
annihilation flux. 

In order to investigate the sensitivity of tidal heating 
to the structure of the micro-halo, we have run additional 
simulations with different initial profiles. We use the cuspy 
profile given by equation 1151 and vary the slope 7 = 0, 



Figure 10. The relative flux of annihilation products from differ- 
ent micro-haloes after a given number of encounters. The typical 
mass loss from a halo would lead to a decrease in flux of between 
a factor of two or three. 



0.5, 1, and 1.5, in addition to our fiducial value, 7 = 1.2. 
This suite of simulations is carried out using a fixed impact 
parameter, b — 0.02 pc, and other parameters as in the 
fiducial run. Our calc ulations are similar to ear lier studies of 
impulsive heating e.g. lAeuilar fc White! <1986ft . who studied 
the structural change in systems with de Vaulouleur density 
profiles. 

Fig. llll shows that up to 30% more mass is lost from the 
cored halo (7 — 0) compared to the cuspy haloes (7 > 1). 
The effect is even stronger for the fraction of mass removed 
from within the scale radius, r s : 2.5 times more material 
is lost from the cored halo. The strongly-bound material 
within the cusp is more stable against tidal disruption than 
that in cored profiles, which have been typically considered 
in previous studies of tidal heating. 

In Fig. |H1 we show the mass loss of the micro-halo for 
ten successive encounters with exactly the same impact pa- 
rameter. This is unlike the situation in our Galaxy where 
micro-haloes orbit the disc for 10 Gyrs near the solar radius. 
It spends about 0.1 Gyr moving through the disc encoun- 
tering stars at a relative velocity of approximately 300 km/s 
and a range of impact parameters. In order to model this be- 
haviour more precisely, we use a Monte-Carlo method, which 
estimates the total amount of mass-loss it would suffer. We 
draw encounter impact parameters from a random distri- 
bution and calculate at each time the stripped mass. After 
each encounter, the density profile of the micro-halo changes 
self-consistently according to the results found earlier with 
the iV-body simulations. 

The mass loss due to the first encounter can easily be 
calculated using equation 1181 From the second encounter 
onwards, the halo density profile changes and it is harder to 
strip during subsequent encounters as we have seen in Fig. 

m 

The mass M a , the halo has after the ath encounter, 
determines the reduction of the mass, which is stripped in 
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Figure 12. Probability density distribution function of the final 
masses of different haloes 



Figure 11. Total mass loss from the micro-halo as a function 
of the inner density slope, 7, determined using the the position 
criterion, r > i? max , after the potential relaxation {triangles), 
as well as the fraction of mass lost from inside the scale radius, 
r < r B (squares). 



the a + 1st encounter. This is unfortunately not a variable 
of equation (1191) . This is a function of n only. So we calcu- 
late at the beginning of each encounter a "virtual encounter 
number" n, which is basically the inverse of equation l|19|l : 



B /log c M a 



.4 



(22) 



For a given M a it computes the corresponding num- 
ber of "standard encounters" with the impact parameter 
b = 0.02pc from figure [5] Now we can calculate the mass 
M a +i after the a + 1st encounter. This must be a deviation 
from equation 1181 with some weighting function w(n,M a ): 



M a 



w(n,M a ) 



1 + 2.1 



(23) 



This weighting function is the fraction of mass, which is 
stripped in the n + 1st standard encounter divided by the 
fraction of mass stripped in the first standard encounter: 



( ' a) "^(D(l-^(n)) 



(24) 



The —^-{n) can be calculated according to equation (1191) . 
This weighting function reproduces the results from Fig. |5] 
Keeping in mind, that [l — ^p(i)] = Ma, equation 123> 
reduces to: 



M a+ i = M a 



M a - exp [A(n + 1) B ] 







-1" 


M(i) 

m y^i 


1 + 2 - 1 (^) 





(25) 



We use this equation recursively for each encounter and then 
repeat the calculation to obtain the probability distribution 
of final masses in Fig. 1121 The central density (at our soft- 
ening length) of a perturbed halo decreases only by a factor 
of about 5, whilst the total mass decreases by an average of 
90 % (see Fig.®. 



5 CONCLUSIONS 

We have studied the disruption of dark matter micro-haloes 
by stars and other substructures using both analytical im- 
pulse approximation and self-consistent iV-body simula- 
tions. The analy tic calc ulations pre sented here are quite 
similar to those of lGreen fc Goodwin! ll200fj) and we come to 
similar conclusions. Our calculations differed in that we used 
more realistic cuspy n-body models and we studied how the 
internal struc ture of these sy s tems evolve due to perturba- 
tions. See also lAngus fc Zhaa d2006T) for a n independent and 
compl ementary study. Earlier studies e.g. lAguilar fc White! 
(1986), also studied cuspy systems and found similar ro- 
bustness to tidal heating in the central regions as we find. 
However our resolution allows us to study the response of 
CDM haloes deep within their central regions. 

• The impulse approximation predicts that those micro- 
haloes in the solar vicinity which formed after z — 130 (about 
80% of the local micro-halo population) should lose most of 
their mass due to close encounters with disc stars. 

• Numerical simulations of individual encounters demon- 
strate that the usual condition of disruptive heating used 
in analytical studies, AE = E^, does not lead to complete 
dissolution of haloes with cuspy density profiles. For the in- 
ner logarithmic slope 7 = 1.2, on average only 30% of the 
mass is lost from the system for this energy change. The 
relation between the fractional mass loss and the energy in- 
put in the tidal approximation is given by equation 1181 : 
AM/M » (1 + 2(A£;/£' b )- 1 )- 1 . 

• The change of particle energies, after the system settles 
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into a new virial equilibrium following the tidal encounter, 
is described accurately by the extension of the impulse ap- 
proximation accounting for virial oscillations. An apparent 
resistance to tidal heating of the material deep in the cusp 
is due to the high binding energy inside the cusp. 

• Repeated tidal encounters lead to diminishing mass loss 
from the same micro-halo. After 10 identical encounters at 
impact parameter b = 2_R v ir, 10% of the halo still remains 
self-bound even though AE = 5E^- 

Near the solar radius within the Galaxy most of the 
mass of the micro-haloes is tidally removed. This material 
forms cold streams in phase space providing a unique signal 
for direct detection experiments. The dense cuspy cores of 
these haloes survive reasonably intact, although the mass 
loss leads to a reduction in annihilation products of about 
a factor of only two to three. These cores could be distin- 
guished by their high proper motions on the sky of the order 
arc-minutes per year. 
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